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A difference of several tenths of a percent has been observed between the direct 
CP asymmetries of D° — > K + K~ and D° — > 7r + 7r~. It has been noted recently 
that CP asymmetries in such singly-Cabibbo-suppressed (SCS) decays can af- 
fect the determination of the weak phase 7 using the Gronau-London-Wyler 
method of comparing rates for B + — > DK + and B~ — > DK~, where D is a 
superposition of D° and D decaying to a CP eigenstate. Using an analysis of 
the CP asymmetries in SCS decays based on a c — > u penguin amplitude with 
standard model weak phase but enhanced by CP-conserving strong interactions, 
we estimate typical shifts in 7 of several degrees and pinpoint measurements 
which would reduce uncertainties to an acceptable level. 
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I Introduction 



The precise determination of phases of the Cabibbo-Kobayashi-Maskawa (CKM) matrix is 
crucial to the understanding of CP violation. At present there appears to be reasonable 
agreement on magnitudes and phases of CKM matrix elements [U[2]. However, discrepan- 
cies among different determinations of these quantities can signal new physics, for example 
due to new heavy particles entering into loop diagrams. Consequently, it is essential to 
pursue the widest variety of measurements of CKM elements. 

One quantity which is determined indirectly with reasonable accuracy but whose direct 
measurement has lagged with respect to many others is the weak phase 7, related to CKM 
elements by 7 = Arg(— V* b V uc [/V* b V c d)- A promising method for measuring 7 directly, 
proposed by Gronau, London, and Wyler (GLW) |3J, compared rates for P> + — > DK + and 
B~ — > DK~, where D is a superposition of D° and D° decaying to a CP eigenstate. In the 
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initial formulation of this method, CP violation in charm decays was assumed negligible, 
as suggested by standard model (SM) estimates [3]. 

Variants of the GLW method include B ± Dif*^ £>*(->■ Dir°, Dj)K ± and B° -> 
DK*° where D decays to CP eigenstates, and processes of this kind in which D° and D° 
decay to a common flavor state such as K~tt + [5 J or to a three-body self-conjugate final 
state such as K s ti + tt~ [5J. Results obtained in these processes have been reported by the 
BaBar Belle jSlHD], CDF HQ, and LHCb [TJl EU El US] collaborations. 

A value of several tenths of a percent has now been seen for AAcp, the difference 
between the CP asymmetries of D° — > K + K~ and D° — > 7i + tt~ [161 02] • These two asym- 
metries have been included in a recent experimental study by LHCb of B — > DK^ |13j . 
concluding that their effect on determining 7 is marginal at the current level of experimen- 
tal precision. Refs. [TSJ [Tj5] have shown that unless CP violation in D decays is taken into 
account, the determination of 7 via the GLW method can be shifted from its true value by 
up to several degrees. 

Three of us have previously assumed that CP violation in charm decays is due to a 
penguin amplitude with the SM phase but enhanced by CP-conserving strong-interaction 
effects [2TJ1 |2"T | |2"2"] . The possibility of such an enhancement was pointed out some time ago, 
in analogy to the likely enhancement of penguin amplitudes in K — > 2n decays [221 [2IJ • A 
large number of authors have suggested that asymmetries at the level observed in AAqp 
cannot be excluded within the CKM framework |25j. A consistent description of SCS 
CP- violating charm decays was found in Refs. [20j [211 122], and predictions were made for 
correlations between CP asymmetries in several decays of charmed mesons to pairs of light 
pseudoscalar mesons P. 

In the present paper we apply this description to B ± — > DK ± ,D — > tt + 7t~ , K + K~ , 
studying within the CKM framework shifts in 7 due to CP violation in D — > PP decays. 
We find shifts of up to a few degrees as noted in Refs. [T8[ IT9] , where no specific scheme 
has been used for CP violation and no correlations existed between CP asymmetries in D 
decays. We identify the most crucial measurements for reducing the uncertainty in effects 
of these shifts to acceptable levels, i.e., below effects of other uncertainties. 

We review the GLW method in Sec. HT1 allowing for CP violation in charm decays, using 
observables involving minimal systematic uncertainties. Present information on direct CP 
asymmetries in SCS D decays to two pseudoscalars is summarized in Sec. IIII1 The approach 
of Ref. [20J is then outlined in Sec. [TV] and applied to obtain predictions for shifts in 7 in 
Sec. [V] We summarize in Sec. ED 

II GLW method in presence of charm CP violation 

We shall be concerned with a single source of information on 7: the decays B ± — > DK ± , 
with D decaying to a CP-even eigenstate such as fo = 7r + 7r~ or K + K~ . Our notation will 
follow that of Ref. [T9]. We define 



A(B~ D K~) 
A(B~ D°K-) 
A(D° f D ) 




(1) 
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where we have taken a strong phase to be zero in B — > D°K with no loss of generality. 
By CP conjugation we then have 

A(B + -)• D°K+) = A B 

A{B + -)• D°K + ) = A B r B e l{5B+ ^ (2) 



A(D^fn) = A f . 

The parameters r B and S B are measurable by combining information from B — )■ DK ± , 
where neutral .D mesons decay to CP-eigenstates, flavor-specific states or Ks^ + ^ ■ Current 
values taken from Ref. [T] are: 

r B = 0.099 ±0.008 , 6 B = (110 ± 15)° . (3) 

The magnitudes of the amplitudes \Af\ and are measurable through the CP- 
averaged branching ratio for D° — > f and £)—>•/ (giving |A/| 2 + |^4/| 2 ) and the direct CP 
asymmetry 

^cp(f) = \ A /!~\i A l • (4) 



We define the weak (CP-violating) phase a/ = Aig(Af/Af). Its value is, as yet, unspecified. 
If <f) is taken as the phase of A f we have 

Af=\A f \e*+; A f = \A f \e^~^ . (5) 

Using the above definitions we may now construct the following amplitudes: 

A(B~ —7- fuK~) = A B A f + A f A B r B e^ &B -^ 

= A B (\A f \ + \A f \ r B e i(S B -j-a f ^ ^ ; 

A(B + ^f D K+) = A B A f + A f A B r B e^ B+ ^ 

= A B (\A f \ + \Af \ r B e l{5B+1+a ^) e^'^ . (7) 



The squared magnitudes of Eqs. (jSJ) and (J7J) give [19] 

\A(B~ f D K~)\ 2 = \A B \ 2 (\A f \ 2 + r 2 B \A f \ 2 

+2 r B \A f \\A f \cos{5 B -l- a f j) , (8) 



\A(B + ^ f D K+)\ 2 = \A B \ 2 (\A f \ 2 + r 2 B \A f \ 2 

+2 r B \A f \\Af\cos(5 B + 7 + "/)) • (9) 



Adding and subtracting the above equations we may form quantities that are relevant in 
constructing the GLW observables: 



\A(B~ ^ f D K~)\ 2 + \A(B + ^ f D K + )\ 2 = \A B \ 2 (\A f \ 2 + \A 



f\ 



l + r| + 2r B cos5 B cos(7 + a / )Vl-(^(/)) 2 j , (10) 
r)| 2 - \A(B+ f D K+)\ 2 = \A B \ 2 (\A f \ 2 + \A f \ 2 ) 
[Af P {\ - r|) + 2 r B sin 8 B sin( 7 + a,) y/l - (i4&(/)) 2 ) • (11) 
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The last expression differs from a similar one in Ref. [19] by a term (1 — r 2 B ) (\Af\ 2 — |v4/| 2 ), 
which vanishes only in the absence of direct CP violation in charm and hence cannot be 
neglected. 

We now take the expressions for Af and Af to be [19] 



A f = \A)\(l + r f e l{& f-^) , A f = \A° f \ (1 + r /e i(<5 ' +7) ) , (12) 

where A® is the amplitude in the absence of a CP- violating term, and we have assumed that 
the source of CP violation has the SM phase —7 as in Ref. [20]. The direct CP asymmetry 
[Eq. 01])] is then given by 

asm) = ; +?Tf M ° 7 ■ (13) 

1 + rj + iTf cos j cos 7 

The fact that Acp(f) is of order a few times 10" 3 in D° ->■ 7r + 7r~ and D° -)■ K + K~ (see 
Tables HI and ITT1 in Section HTB suggests that 77 in these processes is also of this order as 
argued in Refs. [201 125]. Our subsequent analysis (in particular the discussion in Sees. 
|IV] and [V} applies also to larger value of 77, for instance of order 10~ 2 , with an order of 
magnitude further suppression provided by sin Sf. We ignore the fine-tuned solution where 
rf is large, but the strong phases in these decays are of order 10~ 3 . One can show that 

_, / 2rf cos 5f sin 7 + r? sin 27 \ 

af = —tan 



\1 + 2rj cos 5/ cos 7 + r 2 cos 27 J 
w -A^(/)cot5 /; (14) 

where the last approximation holds to leading order in 77. 

It has been suggested in Ref. [26] that when applying the GLW method one normalizes 
the CP- averaged rate for B ± — > }'dK ± by that for B ± — > /d^, 



0/ _ nB-^f D K-) + T(B+^f D K + ) 

«K,,- T{B -^ fD7l - )+T{B+ ^ fDn+) • 



and one takes the ratio of this fraction and a corresponding fraction for D flavor state, 

, , s T(B~ -> D°K~) , , 

fl < A >' s r (B - =; ^-) • (16) 

Significant experimental systematic uncertainties cancel in these fractions [7] El [TT], [T3] . 
Defining ratios of amplitudes and strong phases in B~ — > Dir~ in analogy with B~ — > DK~ , 



^ D 



r BW e l(5s w- 7) , (17) 



A(B~ DH- 
one may express the double fraction 



R K 

Rf ° p+ = r(k^) (18) 



in terms of 7 and these parameters. 
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In the absence of CP violation in D° — > fo, one has [26] 

f l + r| + 2r B cos<5 B cos7 
rLfip, — - 5 , I J-" J 

1 + r z B(n) + 2 r B{7r) cosd B{lT) cos 7 

where the parameter t B {-k) is expected to be very small, r B ( n ) ~ r B ta.ia 2 9c ~ 0.005. [See 
We note that in the approximation of neglecting CP violation in D° — > f B the 
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ratio R s K/n may be defined as R f K/n = [T(B~ D CP+ K~) +T(B + D CP+ K + )]/[T(B- -4 

+ — > D C p+7T + )]- Consequently this ratio and the double ratio R C p+ do 
not depend on f D . 

Including CP violation in D° — > one finds an expression for Rq P+ which depends on 
fo through the CP asymmetry A^ P (f), 



l + rl + 2r B cos5 B cos( 7 + a f )y/l - {Af P {f)) 2 
Kq P , = , . (20 ) 

1 + r| (w) + 2 r B(n) cos5 B{w) cos( 7 + a f )y/l - (A&(/)) 2 

Neglecting corrections in R CP + which are quadratic in A^ P (f) [C(10~ 5 )] and using Eq. (JHJ), 
we note that corrections linear in A^f p (f) ~ fewx 10~ 3 are multiplied by r B and are therefore 
negligible relative to R CP + = 1 + 0(r B ). Thus, a comparison of Eqs. ( !T9|) and (|20|) shows 
that the determination of 7 is not affected in a significant way by including A^ P {f) in 



r! 



cp+- 

The other measurable quantity used in the GLW method is the CP asymmetry A(, p+ , 

A f _ TPS' -± fpK~) - T(B+ -» f D K+) 
Acp+ ~ T{B- -> fuK~) + -> • 

In the absence of CP violation in D° — )■ this asymmetry may be defined as [r(l? _ — >■ 
Dcp+K-) - T(B + ->■ J D CP+ K+)]/[r(5- Dcp+K-) + r(5+ ->■ D C p+^ + )] which is 
independent of /d and is given by [26] 



j 2 rg sin 5 B sin7 

1 + r B + 2 rs cos cos 7 



(22) 



When including CP nonconservation in D° — > fp> the asymmetry Aq P+ becomes dependent 
on Ac P (f). Neglecting terms quadratic in A(fp(f), one finds 

A f _ Af P {f){l -r 2 B ) + 2 r B sin S B sin( 7 + a f ) 
1 + + 2 r B cos os cos(7 + «/) 

The two terms in the numerator and the last term in the denominator involve corrections 
linear in A^ P (f) modifying the expression (122"]) for the case of no direct asymmetry. 

We note that an expression independent of A^ P (f) similar to (1221) . but with an op- 
posite overall sign and an 0(r B ) correction with opposite sign, describes the asymmetry 
for Cabibbo-favored D° decays to CP-odd eigenstates (such as Ks4> and Ksir ) where CP 
violation is negligible in the CKM framework, 



f 2 r B sin 5 B sin 7 

1 + r B — 2 r B cos B cos 7 



(24) 
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A precise measurement of this asymmetry could avoid uncertainties from A^ P (f). 
Writing 

2 r B sin 5 B sin j eS 



A f 



1 + r B + 2 r B COS 5 B COS 7eff 



(25) 



we assume a measurement of Aq P+ from which 7 er r is determined. Our purpose is then to 
study the shift £7 = 7 — 7 er r as a function of A^ P {f). Defining 7 + a/ = 7 or T + £7 + a/ = 
7efr + w e expand the numerator and denominator in (123]) to first order in x. Comparing 
this expression with ( 125]) . cross-multiplying the two ratios, cancelling the leading terms and 
keeping terms linear in x and AQ P (f), one finds 

1 + r% + 2r B cos 5 B cos 7 er T 



1-ri 



A dir 



(/) 



cot 5 4 



2r B sin5 B (1 + r|) cos7 cff + 2r B cosS B 
1 + r% + 2r B cos 5 B cos 7 eff 



2rg sin 5 B (1 + r|) cos 7 cr r + 2r B cos <5 B 



2r/ cos5/sin7e ff - A$p(/) 



2r B cos 5s cos7eff 



(26) 



_2r B sm5 B (1 + r|) cos7 eff + 2r B cos5 B _ 
While both corrections are exhibited as proportional to A Q P (f), the first term appears to 



diverge as 5f — > 0. Recalling that this correction is actually ~ 2r/cos(5/sin7 e ff [Eq. (113])]. 
which is finite at 5f = 0, we will use the last expression in Sec. |V] The second term in ( 126]) 
illustrates an enhancement by l/2r B of the shift in 7 due to A^ P {f) [T8] [T9]. 

We mention in passing two approximations which have been applied to the CP asym- 
metry in (123]) . Neglecting 07 [which is of the same order in 77 as A^p^f), but leads to 
a correction linear in r B ] and approximating the overall coefficient of A^ p (f) by one, Eq. 
reduces to an expression employed in Ref. 



A f 



2r B sin S B sin 7 



+ ^cpU) ■ 



1 + r B + 2r B cos S B cos 7 



(27) 



(28) 



A further approximation, keeping only 0(r B ) in the first term, leads to 

A f CP+ ~ 2r B sin 5 B sin 7 + A A j P (f) . 

A CP asymmetry in B ± /dK^ has been measured in Refs. [7] [TT] [12] consistent with 
an estimate Aq P+ (B ± — > f^K^) ~ 2r# sin #5 sin7 ~ 0.15 — 0.20. The current error in the 
world-averaged value [28], Acp+ = 0.19 ± 0.03, is still too large for sensitivity to A^p^f). 
The corresponding error in the world-averaged measurement, Acp- = — 0.11±0.05 favoring 
an opposite sign as anticipated, is somewhat larger. 

A much smaller asymmetry is expected in B ± — > Dtt^ where we estimated r B M ~ 0.005. 
This asymmetry is given by an expression similar to 



A f CP+ {B -> f D n) 



A cp(f)0- ~ r B(n)) + 2 r B(n) sinS BM sin(7 + aj 

1 + r l(7r) + 2 r B(7r) cos 8 B (tt) cos(7 + a f ) 
2 r B(7T) sm5 B{7v) sin 7 + A d J P (f) . 



(29) 



The two contributions, which in principle may be disentangled by measuring also Acp-{B — > 
/ D 7r) ps —2 rs(7r) sin5s(7r) sin 7, are of comparable magnitudes, each less than a percent. 
Thus, while the CP-averaged rate for B ± — > is suitable for normalization [see Eq. 

(II 5p ], a useful measurement of the corresponding asymmetry does not seem feasible in the 
foreseeable future. 
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Ill Present information on Afj T P 

Dedicated measurements of the difference AA CP = A^p^D -)• K + K~) - A^pi^D -)> 
7T + 7r _ ) (in which many systematic errors cancel) have been performed by the LHCb [16J 
and CDF [17] Collaborations, while Belle has combined independent measurements of the 
two asymmetries [27] to obtain a value of AAcp- The results are shown in Table [H We shall 
assume that measured asymmetries are equal to direct ones, neglecting possible contribu- 
tions from indirect (mixing-induced) asymmetries which would lead to slightly different 
averages [H [28] . 

Table I: Experimental results on AA CP = A^ P (D° -»■ K + K~) - A$ P (D° tt+tt-). 



Reference 


Value (%) 


LHCb [TB] 


-0.82 ±0.21 ±0.11 


CDF [T7J 


-0.62 ±0.21 ±0.10 


Belle [27] 


-0.87 ±0.41 ±0.06 


Average 


-0.74 ±0.15 



Table II: Experimental results on some direct CP asymmetries in D decays. 



Decay 


Reference 


Value (%) 


D° 7T + 7T- 


CDF [29] 


0.22 ±0.24 ±0.11 




Belle [27] 


0.55 ±0.36 ±0.09 




Average 


0.33 ±0.22 


D° -»■ K + K- 


CDF [29] 


-0.24 ±0.22 ±0.09 




Belle [27J 


-0.32 ±0.21 ±0.09 




Average 


-0.28 ±0.16 


D + -> K + K° (a) 


BaBar [30] 


0.46 ±0.36 ±0.25 




Belle [27J [31] 


0.08 ±0.28 ±0.14 




Average 


0.21 ±0.25 



(a) After subtraction of CP asymmetry due to K°-K mixing. 



The average in Table [I] will be used in the next Section to constrain the magnitude of a 
SM penguin amplitude as a function of its strong phase. (Slightly different averages were 
used in Refs. [?] and [19].) In addition weak constraints on this strong phase will be seen 
to result from measured CP asymmetries in individual final states, quoted in Table HTl 

IV Charm CP violation with enhanced SM penguin 

In Ref. |20j three of us have calculated CP asymmetries for several D — > PP decays, where 
P = 7r, K, assuming that the nonzero value of Aq P is due to a SM penguin amplitude with 
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Table III: Magnitudes Tf and phases (j) T for some D — > PP processes. 



1 loprnr 

uecay 




<Pt — Ar Sv J /J 


mode 


(lO" 7 GeV) 


(degrees) 


D° tt+tt- 


4.70 


-158.5 


D° K + K~ 


8.48 


32.5 


D° -> 7T 7r° 


3.51 


60.0 


D + -> TT+F 


6.87 


-4.2 



the weak phase of the standard model c — > b — > u loop diagram, but with a CP-conserving 
enhancement as if due to the strong interactions. In this case the magnitude and strong 
phase of this amplitude Pb are correlated in order to fit the observed CP asymmetry, 
allowing the prediction of CP asymmetries for other singly-Cabibbo-suppressed modes. 
We refer the reader to that work for details, but outline the method briefly here. 

For the decay of a charmed meson D to any final state / we are defining a direct CP 
asymmetry using the same convention as in (J4]) 

^ m _ rp -> /) - r(5 -> /) 

We take the CP-conserving amplitudes from a flavor-SU(3) description of charm decays 
presented previously [321 [33] . A new ingredient [20] with respect to that work is a U-spin 
breaking c — > u penguin amplitude whose V* d ^V u d( s ) term reproduces satisfactorily decay 
rates for singly-Cabibbo-suppressed (SCS) processes, while its V* h V u b term accounts for 
AAqp- This scheme allows the prediction of other CP asymmetries in SCS charmed meson 
decays. 

Denoting by 5f a phase defined in Ref. J2U] as <f>f, we write the amplitude for a decay 
D — > f in a manner similar to (I12p 

A(D -)>/) = A f = |T>|e^(l + r /e i{ ^" 7) ) . (31) 

Here Tf represents terms with the weak phase of the tree-level terms contributing to that 
amplitude, <pj, is its strong phase, Tf is the ratio of the magnitude of the CP- violating 
penguin contribution to that of Tf, —7 is the weak phase of the CP- violating penguin, and 
Sf is the strong phase of the CP- violating penguin relative to Tf. Here one has 

5/ = Arg(P 6 )-4 + 7 , P h = pe ^\ (32) 

leading to the relation 5f = 5 — <^ T . The magnitudes Tf and phases (fij. for some D — > 
PP processes [20] are summarized in Table IIHI For D° — > 7r°7r°, the amplitude must be 
multiplied by an additional factor of — 1/\/2 [20J. 
The CP asymmetries are then 

Af P (f) - 2r/ sin 7 sin 5/ 



1 + rj + 2rf cos 7 cos Sf 

2p\Tf\ sin7 sin(5 — 0j.) 
\T f \ 2 +p 2 + 2p\T f \ cos7cos(<5 — 



(33) 
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s 

Figure 1: Magnitude p of the CP-violating penguin amplitude as a function of its strong 
phase 5. The red curve was obtained using the value AAcp = —0.74%, while the inner 
(blue) and outer (green) bands respectively correspond to ±lcr and ±1.64(7 shifts from this 
value, where a(AAcp) = 0.15%. The plot is shown only for 7 = 67.2°, as p is not very 
sensitive to the exact value of 7. 

Taking 7 = (67.2^^)° from Ref. pQ, the world-averaged asymmetry AAcp = (—0.74 ± 
0.15)% from Tabled] is used to constrain the magnitude p of the penguin amplitude as a 
function of its strong phase 5. The value of p, plotted in Fig. [TJ is nearly constant at several 
tenths of a percent of the amplitudes in Table ITTTl for a wide range of 5. This corresponds to 
values of r f of this order. The value of p and corresponding values of are about an order 
of magnitude larger for extreme values of 6 in Fig. [TJ The constraint on p as a function of 6 
allows one to predict asymmetries for (e.g.) D° — > 7r + 7r~, K + K~, tt°tt° and D + — > K + K°, 
as plotted in Fig. [2] for 7 = 67.2°. Very similar results (not shown) are found for 7 = 71.6° 
and 62.6°. 

For much of the range of 5, Aq P (tt + tt^) is predicted to be positive while Aq P (K + K~) 
is predicted to be negative. (In the U-spin limit they would be equal and opposite.) This 
is consistent with the central values in Table [TTJ However, the predicted central values 
for A'^p(K + T{ ) are negative for much of the range of 5, whereas the 2a lower limit 
A'^p(K + l( ) > —0.3 would tend to favor values of 5 > —tt/2. Improved measurements 
of all these individual CP asymmetries would of course be highly desirable. 

V Predictions for shifts in weak phase 7 

Taking the predicted values of Aq P , Eq. (|26|) implies a shift in the weak phase 7 associated 
with the use of each D-decay process with a CP-even final state. In Fig. Owe present these 
shifts for the final states n + 7r~ and K + K~ . In addition to the central value of shifts £7 we 
also present errors in £7 due to the variation of the various measurable parameters. In our 
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Figure 2: Direct CP asymmetries for some SCS D — > PP' decays. Curves and bands as in 
Fig. HI 7 = 67.2° is assumed. Very similar plots (not shown) are obtained for 7 = 71.6° 
and 62.6°. The dashed horizontal (red) lines in the lower right panel denote 90% C. L. 
limits based on the average of Belle [271 E] an d BaBar [30] results. 
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Figure 3: Shifts £7 as calculated from Eq. (126]) for D° — > 7t + tt~ (left panel) and D° — > 
K + K~ (right panel). The solid central (black) curves denote the central value for £7. Also 
shown are ±la errors in rp (short dashes [blue]), 5b (medium dashes [blue]), and AAcp 
(long dashes [blue]). The (red) dots represent the effect of ±lcx shift in the measured value 
of 7 on the central (black) curve. The solid outer (red) curves denote the effect of a ±lcr 
shift in the measured value of 7 on the curve that is obtained by adding in quadrature the 
errors in re, 5b, and AAcp- Here 7 = (67. 2^4^)°. 

calculations of the shifts £7 we have used the following values of r#, 5b, and 7 taken from 
Ref. [I]: 

r B = 0.099 ± 0.008 , 8 B = (110 ± 15)° , 7 = (67.2t££)° . (34) 

In Fig. [3] the central black curves represent £7 as obtained from Eq. (|26|) using central 
values of the measurable parameters rg, 5b, AAcp an d 7- The errors in £7 from the 
la error in the measurement of the first three parameters are shown in blue using short, 
medium, and long dashes respectively. In order to obtain the effect of varying 7 on the 
shifts £7 we use Eq. (T26]) with the value of 7 set to its drier limits given in Eq. (134"1) . while the 
other parameters are held fixed at their respective central values. This effect is represented 
by the red dots in Fig. [31 In order obtain the overall error in £7 we first add in quadrature 
the errors due to r B , 5b, and 7. We then estimate the effect of varying 7 between its ±la 
limits on the quadrature sum. This effect is represented by the solid red curves in Fig. [3J 

We see in Fig. [3] that the central value of the shift in 7 for the ir + n~ state is nearly 
constant over the entire range of allowed values of the strong phase 5. This appears to be 
the effect of an accidental cancellation between a variation of A^ r p which is modest to begin 
with and a compensating factor due to the variation in 5f for 7 = 67.2°. When the various 
sources of errors on #7 are included, this no longer holds, as depicted by the la boundary red 
curves in Fig. [3j On the other hand, £7 obtained using the K + K~ state shows appreciable 
variation (~ 7°) over the allowed range of 5. A similar exercise when performed using 
the 7r°7r° state yields an even larger variation in £7. However, rate asymmetries involving 
multiple neutral pions in the final state are not expected to be measured with adequate 
precision in the foreseeable future. We have therefore chosen to omit the 7r°7r° decay mode 
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from the present discussion. 

In view of the variety of dependences of direct CP asymmetries on 5, and different 
behavior of shifts in 7 for 7r + 7r~ and K + K~ final states, it would be beneficial to further 
pin down 5 [e.g., with a better measurement of A^ P (K + K )], and to apply the GLW 
determination of 7 to the widest assortment of CP eigenstates in neutral D meson decays. 
After the appropriate shifts £7 have been taken into account, inconsistencies in final values 
of 7 obtained for different charm final states could point the way to effects of new physics. 

VI Discussion and summary 

The determination of the weak phase 7 by means of the CP asymmetry in DK ± , 
followed by the decay of the neutral D meson to a CP eigenstate [3], must take account 
of asymmetries in the decays of the neutral D mesons [131 HI EH]- We have calculated 
the corresponding shifts £7 in an approach which imagines these CP asymmetries as due 
to a c — y u penguin amplitude with weak phase of the standard model but enhanced by 
(presumably nonperturbative) strong interaction effects beyond those anticipated by the 
majority of authors. The observed value AA C p = A$ P (D° -)> K + K~) - A$f P (D° -)■ 
7r + 7r~) = (—0.74 ± 0.15)% has been taken as a constraint, leading to a correlation between 
the magnitude p of the CP-violating penguin and its strong phase 5. 

For 7 = 67.2° the central value of the shift associated with the decay D° — > 7r + 7r~ 
is approximately 5 , j(tt + 7i^) ~ —4.4°, with little dependence on the allowed range of the 
strong phase 5. However, the factor multiplying A^p(/)/2r_B in the last expression for £7 
in Eq. ( 126|) is roughly inversely proprtional to cos 7, which is fairly small and fairly sensitive 
to 7. Thus, when 7 is varied within its currently allowed range of about ±4.5°, the value 
of 57 varies considerably. It is further affected by uncertainties in AAcp, re, and 5p, and 
acquires some dependence on 5. 

The shift associated with D° — > K + K~ is of the other sign (as is Aq P ) and depends on 
both AAcp and 5 as well as the uncertainties in r# and 5b- Using measurements of Aq p 
for both these two decays, with the help of improved knowledge of A^f P {D + — > K + K°) to 
pin down 5, the uncertainty in 7 due to CP violation in charm decay can be reduced to a 
level where it is no longer the dominant uncertainty when applying the GLW method to 
the decays B ± ->■ DK ± . 

Let us be clear about the limitation of our study. The shift in 7 we calculate using 
our c — y u penguin amplitude model is based on measuring Rcp+ and Aqp+ in B~ — > 
(it + -n~)p)K~ , B~ —7- (K + K~)p>K~ and their CP conjugates, taking rp and 5b as given. 
This assumes that one has measured r# and 5b first and then uses only these GLW processes 
to determine 7. 

An actual analysis for determining r B , 5b, and 7 in B T — > DK T [HI [TOl [TTJ [T2j [13] 
combines information from D decays to CP eigenstates [3], flavor states [5] and three- 
body self-conjugate final states [Bj. In this global analysis 7 is also constrained by rates 
and asymmetries in B — > DK where there is no direct CP violation in D decay, for 
instance in decays to flavor states and CP-odd eigenstates. [See Eq. (I2"4"j) ]. Thus any 
actual determination of 7 from B — > DK will involve a considerably smaller shift than we 
calculate. 

Furthermore, it has been noted [26l 1341 [35] that in the self-tagged decays B° — > DK*° 
the ratio r* B of B° -»■ D°K*° and B° -»■ D°K*° amplitudes, both of which are color- 
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suppressed, is expected to be about three times larger than tb defined in B + — > DK + . 
We have seen an enhancement by l/2rg of the shift in 7 due to A^p(7i + 7i^ , K + K~) in 
B + — > DK + . This implies that when applying the GLW method to B° — > DK*° the shift 
in 7 due to these direct CP asymmetries in D° decays is expected to be about three times 
smaller than calculated above. First measurements of relevant observables in B° — > DK*° 
and its CP-conjugate have been reported very recently by the LHCb collaboration [T5] . 
Early measurements of these processes have been performed by B aB ar [36J . 
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